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1. Introduction and methods 

Recent years have shown remarkable progress in lattice simulations of quantum chromody- 
namics (QCD). Dynamic simulations with multiple flavors of light quarks are now routinely per- 
formed and there are first calculations at physical pion masses. Properties of ground states can 
now, in many cases, be determined with small statistical errors and with fully controlled system- 
atic uncertainties. At the same time new computational methods for smeared interpolating fields 
[jl|, and the use of the variational method [|]-|5]] enable extraction of excited energy levels both 
in an unprecedented number and with unprecedented statistical accuracy However, most 

hadronic excitations are resonances and can decay strongly. Therefore the interpretation of the re- 
sulting data is straightforward only in the limit of very narrow states, which is often not the case in 
nature. 
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T [MeV] 


6i(1235) 
fli(1260) 


142 ±9 
250 - 600 


£*(1410) 
*S(1430) 


232 ±21 
270 ±80 


D^(2400) 
Di(2430) 


267 ±40 
384±J30 



Table 1: Examples of light, strange-light and charm-light meson resonances. Values from the PDG compi- 
lation [§]. 



Indeed, taking a look at the meson tables in the Particle Data Group (PDG) compilation [Q], 
resonances with a substantial hadronic width are more the norm than the exception. Table [j] lists 
examples for light, strange-light and charm-light mesons and their PDG values for the resonance 
parameters. As can be seen, there are many examples of resonances with sizable total widths, 
which should not be neglected. In these proceedings, progress in extracting resonance properties 
from lattice simulations is reviewed with an emphasis on recent simulations data rather than the 
associated theory. 

In experiment, resonance properties are extracted using partial wave analysis. For the relatively 
simple case of elastic scattering, the scattering amplitudes a\ is related to the scattering phase shift 
5/ for the 1-th partial wave: 

e 2i5/ _ i 

fl / = sin5 / e , * = — — - . (1.1) 
2i 

Near a single relativistic Breit-Wigner shaped resonance, one can then parametrize the scattering 
amplitude a\ in terms of a resonance position sr = m\ and decay width F 

_ -y/sT( S ) 

S — SR + ly/sF{s) 

In a lattice QCD calculation in Euclidean space there is no direct access to the scattering 
amplitudes [|l0|]. However, as has been pointed out by Liischer [JTT|-[l3[], the phase shift of the 



continuum scattering amplitude in the elastic region can be determined from the discrete spectrum 
in a finite box. 
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L » few fm 




L= few fm 



Figure 1: Illustration of hadrons with momenta p\ and p2 in a finite volume. For a description please refer 
to the text. Figure taken from Jl4| ] . 



1.1 The Luscher method 

As mentioned, there is a relation between the phase shift of the continuum scattering amplitude 
in the elastic region and the discrete spectrum in finite volume JTT|-[T3|]. Figure [j] illustrates the 



idea. In a large box with spatial volume L 3 and L » 1 fm (depicted on the left-hand side) the energy 
of the two particle system is to a good approximation given by 

E = E(p 1 )+E(p 2 ), (1-3) 

where E(p) = \fijn 1 + p 2 ) and p = n(2n/L) in a relativistic simulation with periodic boundary 
conditions in space. In a small box (depicted on the right-hand side) with L ~ 2 ... 5 fm the en- 
ergy of the interacting system is noticeably shifted with regard to the non-interacting energy level. 
This energy shift is related by Liischer's formula to the elastic scattering phase-shift. Extracting 
resonance parameters from a lattice simulation using this relation therefore involves the following 
basic steps: 

(1) Extract the energy levels E n (L) in a finite box for one or more box sizes L. 

(2) The Luscher formula relates this spectrum to the phase shift of the continuum scattering 
amplitude. 

(3) Given the relevant phase shift data one can extract the resonance mass Mr and the width F R 
or the coupling g with some degree of modeling/approximation. 

To illustrate the finite volume spectrum, the few lowest energy levels for a p -meson-like sys- 
tem are plotted in Figure ||. For a small coupling the typical avoided level crossing behavior is 
observed (left panel). With an increased coupling g the energy shifts with respect to the non- 
interacting energy levels get more pronounced and the avoided level crossing gets more and more 
washed out (middle panel). The effect of increasing the mass at fixed coupling is also illustrated 
(right panel). From these illustrations it is clear that data for just a couple of energy levels in a 
single volume is of limited use. For a given lattice ensemble more data can be obtained by also 
considering moving frames. 



While the original relation [jn]-|T3p was limited to rest- frame calculation in multiple spatial 
volumes L 3 , the corresponding Luscher formulae and interpolator constructions for equal mass 
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L nix % 




Figure 2: Illustration of levels in a finite volume L 3 x T for a p-like resonance. The resonance mass m 
in units of the pion mass and the coupling g are varied (assuming only elastic scattering and using the 
appropriate Liischer formula). The blue levels show the non-interacting levels and the blue dotted line 
indicates the resonance mass. 



hadrons mux = mui in moving frames [|15|-|T8|] have been known for a while. During the last year 
the corresponding expressions for unequal mass hadrons n%h\ ^ nihi in moving frames have been 
derived [[jj]-^2|]. It turns out that for the case of unequal mass hadrons, even and odd partial waves 
mix, creating an additional complication for the ab initio determination of phase shifts 5/ from 
lattice data. 

In addition to the above, there are some more recent ideas to extract resonances from lattice 



simulations [23, 24]. One of these, the so-called histogram method [23] is illustrated briefly in the 



next section. For a description of the so-called correlator method please refer to Q24| ] directly. 

2. Methods in a toy model study 

In a recent publication, Giudice, McManus and Peardon [[25|] tested both the Liischer method 
and the histogram method [E3] in the context of the 0(4) non-linear sigma model, where precise 



data can be obtained: the authors of 0250 managed to accurately extract six energy levels for a 
range of L/a = 8, 9, . . . , 19. For both methods, data was extracted in the elastic region (where the 
methods should be applicable) and in the inelastic region. 



The histogram method [ |23| ] can be summarized by the following procedure: 

(1) Determine the lowest few energy levels and interpolate to obtain a continuous function E n (L) 
in the interval [Lo,Lm]- 

(2) Slice this interval [L ,L M ] into M parts of length AL = 

(3) Slice the energy interval [E mini E max ] into bins of length A E . 

(4) Make a histogram and normalize to get the distribution W(E) or correspondingly W(p). 

(5) Subtract the non-interacting background Wq(E) (or Wq (/?)). 



(6) In [23] it is shown that close to a resonance one gets a Breit-Wigner shape 



1 dE 

^'-»»^ BFW ' n P )=nE)j- p . <2.D 
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Figure 3: Probability distribution W — Wo for a narrow resonance with different strategies for the background 
subtraction. Figures taken from [E5n. 



£ o 




Figure 4: Final results for a resonance in the inelastic regime. In this case, the relation to resonance param- 
eters is unclear. Figure taken from [S] . 



For further details please refer to [|23|, |25p. In practice, it turns out that the successful applica- 
tion of this procedure hinges on details of how to deal with the necessary background subtraction 
(step (5)). This behavior is illustrated in Figure ||, where a narrow resonance was investigated with 
different strategies for the background subtraction. The figure also illustrates that, in the context 
of this toy model, a narrow resonance can be extracted reliably once the background has been sub- 
tracted correctly. Unfortunately this was only achieved after omitting some of the data from the 
analysis. 

The authors also apply the histogram method in the inelastic case, where there is no theoret- 
ical support for its applicability. Figure || shows that nevertheless some sort of resonance shape 
emerges. 



Let us now turn to the Luscher method [Jllj-|13p described in the previous section. The left and 
center panels in Figure [5] show the results for a narrow and wide resonance. In both cases a clear 
resonance shape can be identified, although the resonance parameters can be determined much 
more accurately for a narrow resonance. In addition the authors also investigated what happens 
in the toy model when the method is applied above inelastic threshold. This can be seen in the 
right panel of Figure where a narrow resonance is seen below threshold and the data becomes 
nonsensical above inelastic threshold. This illustrates that care needs to be taken when using the 
Luscher method for the extraction of resonance parameters. 

To compare how well the two methods fare, Table lists the resonance parameters obtained 
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Figure 5: Phase shift data obtained from the Luscher method.The blue line in the right plot shows the 
position of the inelastic threshold. For a description please refer to the text. Figures taken from | 





Luscher 


histoj 


*ram 


parameter set 


aM G aT a 


aM 


aT G 


set A 


1.35(2) 0.115(8) 


1.33(5) 


0.10(5) 


setB 


2.03(2) 0.35(2) 


2.01(2) 


0.35(10) 


set C 


3.1(7) 1.2(5) 







Table 2: Comparison of resonance mass M and width T obtained with the Luscher and histogram methods 
for narrow (set A), medium (set B) and broad (set C) resonances [fZ5|]. 



with the help of the two methods. While the results from the two methods agree within their 
statistical uncertainty, the Luscher method leads to smaller statistical uncertainties and can be used 
to successfully extract broader resonances. 

3. QCD resonances 

Let us now turn our attention to the calculation of hadron resonances in QCD. These calcu- 
lations are computationally demanding and very challenging from a technical point of view. It is 
observed in many calculations that qq operators conventionally used for the study of hadron ex- 
citations couple very weakly to multi-hadron states [^6] -^8], §]. A similar observation was made 
in string breaking studies |3(J. This necessitates the inclusion of hadron-hadron interpolators. 
Figure |6] shows an example of quark diagrams for a typical meson-meson channel. In particular 
some of the quark lines contain backtracking quark lines which are very expensive to calculate as 
they require the use of all-to-all propagator techniques. A particularly promising technique is the 
distillation method [|, g]. 

3.1 The p meson: A benchmark calculation 

In experiment, the p meson is seen as a p-wave resonance in 7r;r-scattering. In many ways the 
p meson is an ideal candidate for a benchmark calculation. While the physical p can also decay 
into four pions, there is a large span of (unphysical) pion masses where the decay into two pions 
is the only possible decay. Furthermore, the resonance is of medium width (much smaller than 
its mass) and well isolated from other resonances. The signal in this channel is of a good quality 
even with a moderate number of gauge configurations. Naturally this was the first QCD resonance 
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Figure 6: Example of the Wick contractions arising in a meson-meson channel with qq and meson-meson 
interpolator basis. 




where Liischer's method has been applied to extract the mass and width of the state plQ. In the 
meantime several groups have successfully demonstrated the feasibility of determining resonance 
parameters of the p meson in a QCD calculation []32|-|37|]. In this section the current state of affairs 
is illustrated. 

Figure ^ shows the results obtained by the ETM Collaboration [|32|]. They use four ensembles 
with the same lattice spacing and different light valence quarks and three momentum frames. This 
enables them to extract phase shifts for each ensemble and to perform an extrapolation to the chiral 
limit (assuming the coupling g pK7l to be mass independent), while all of their data is extracted in 
the elastic regime where the framework is applicable. 

Figure || shows similar results from the PACS-CS collaboration using two ensembles and three 
momentum frames [^3]]. Again, all data has been obtained in the elastic region. 

Figure |9] shows results by Lang et al. [j^ from a single ensemble, demonstrating that a very 
good statistical accuracy can be achieved using the distillation technique [[]], §]. All these results 
were presented in parallel talks at Lattice 2011. 

More recently, Pelissier and Alexandru j35| ] presented results from a dynamical simulation 
on asymmetric lattices []38], [39|], which followed an exploratory quenched simulation [[tt]]. They 
generated three ensembles with Nf = 2 nHYP-smeared clover fermions at a lattice spacing of 
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a = 0.1225(7)fm and with a pion mass of m n = 304(2)MeV. Their lattice sizes are 24 2 x L z x 48 
with L z = 24, 32,48. To extract the two lowest states they use a variational basis with one qq and 
one meson-meson interpolating field. These results are plotted in Figure [U]. 

To compare all results for the resonance mass Mr and the coupling g pKK it is best to plot the 
results in dimensionless units by multiplying the masses with the Sommer scale [ pT[ ] r$. To this 
end the values for ro/a for each ensemble are used. The results are shown in Figure |TT] and the 
errors plotted encompass the statistical uncertainty and the uncertainty in the Sommer scale. Plotted 
in this scale-independent way all available results agree with each other fairly well, although the 
values obtained by Lang et al. for both the mass and the coupling are somewhat lower than the 
rest. This highlights the need to address various systematic uncertainties in more detail in future 
simulations. 

In this regard, it is instructive to take a look at the interpolator dependence of the results from 
[Q which are shown in Figure [l2|. On the left side of the figure the energy levels resulting from 
several choices including both pion-pion and qq interpolating fields are shown. Here, a basis of 
just two interpolating fields turns out to give rather unreliable results when compared to several 
choices of a bigger basis. Notice that all fits had an acceptable % 2 /d.o.f.. On the right side of the 
plot choices of basis without a pion-pion interpolator are shown. In all cases the errors are larger 
than for the mixed basis, especially for the first excitation. In the case of a moving frame with 
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Figure 11: Comparison of recent lattice results for the p meson resonance mass (right panel) and coupling 
(left panel). 



P = (1,0,0) the results are even inconsistent for a small basis and only become consistent with 
the mixed basis when the basis is enlarged. It should therefore be stressed that a suitable basis is 
crucial to extract physical results. 

3.2 Recent results for other QCD resonances 

While the extraction of the p resonance is a great proof of principle that QCD calculations of 
resonance properties are feasible, this can only be the start. In this section recent results in other 
channels are reviewed. 

3.2.1 Meson-meson scattering in the Kk, Dk and D*n channels 

Recently, first steps towards the determination of phase shifts and an extraction of resonance 
parameters have been taken for the case of meson-meson scattering in the Kit, Dn and D*7i chan- 

N[xN T n p a[fm] L[fm] #configs mJMeV] m^[MeV] 
16 3 x 32 0.1283 7.1 0.1239(13) L98 280/279 266(3)(3) 552(2)(6) 

Table 3: Parameters for the nHYP-smeared Wilson-clover lattices Q used in Q 
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Figure 12: Interpolator set dependence for the energy levels used to extract the p resonance in Q]. 
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Figure 13: Phase shift results for ^scattering in s-wave and p-wave for isospin \ and |. 



nels []44|, |45|]. For both of these studies configurations with Nf = 2 flavors of nHYP smeared 
Wilson-clover quarks Q were used. Table |5| shows the relevant parameters. The distillation 
method ^ ||] was used to construct an interpolator basis of several qq and meson-meson interpola- 
tors. In both cases only the frame with total momentum zero was considered, which does not suffer 



from mixing of even and odd partial waves [19-22]. 



Figure O shows phase shift results from Kit scattering for isospin A and §, in s-wave and 



p-wave. To illustrate the qualitative agreement, data extracted from experiment ^] is shown 
in black, blue and green while the lattice data is shown in red. For more information and some 
cautionary remarks please refer to [p4|]. 

Figure |l4| shows recent lattice results [^8|-^, for the s-wave scattering lenghts in units of 



the reduced mass /i^ 1 compared to results from Chiral Perturbation Theory (#PT) []5T|] or from a 



Roy-Steiner analysis p52p. 

A similar approach has been used to investigate Dn and D*7T scattering [[l5|] where the charm 
quarks are treated with the Fermilab method [|53|]. Figure [l5| shows ^cot<5/ for the Dq and D\ 



In leading order ^PT this ratio does not depend on m n . 
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Figure 14: Compilation of results for the isospin ^ and | s-wave Kn scattering lenghts taken from [|14|. For 
more detail please refer to the text. 
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Figure 15: Phase shift results for the Dn and D*n channels. To extract the data for D*k scattering the heavy 
quark limit has to be assumed. 



channels, since the combination ^ cot 5/ = p- (s — mj:) is linear in case of a Breit- Winger resonance. 
The left panel shows results in the J p = + channel where the Dq(2400) resonance is observed in 
experiment. Three levels have been extracted. 

In the J p — 1 + channel there are two resonances D\ (2420) and D\ (2430) and without further 
assumptions the data is insufficient to extract resonance parameters. In the heavy quark limit, the 
narrow D\ (2420) is expected to exclusively decay in d-wave, while the broad D\ (2430) will decay 
exclusively in s-wave. In this limit, one then obtains the results displayed in the right panel of 
Figure |l5|, where the level associated with the D\ (2420) has been omitted. 





D£(2400) 


Di(2430) 


g lat [GeV] 


2.55 ±0.21 


2.01 ±0.15 


g exp [GeV] 


< 1.92±0.14 


< 2.50 ±0.40 



Table 4: Resonance parameters for the Z)q(2400) and for the D\ (2430) extracted in Experiment values 
are quoted as upper limits to reflect that only the total width is known. 



Assuming a Breit- Wigner shape for both the Dq(2400) and for the Di(2430), the resonance 
parameters displayed in Table || are obtained. In addition to the channels treated as resonances, 
energy levels were also extracted in other channels. The resulting spectrum is shown in Figure [16. 
For further information including a detailed description of assumptions please refer to 
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Figure 16: Energy differences AE = E — | (M^ + 3M#* ) for Z) meson states on the lattice and in experiment. 
Likely quantum numbers have been assigned to some of the states seen by BaBar [Q]. For more information 
please refer to the text and to [p3[] . 



3.2.2 k and a resonances from staggered simulations 

There are two recent papers [|55|, |56|] aimed at the extraction of the k and a resonances using 
staggered quarks. In both cases the asqtad medium-coarse ensemble of size 16 3 x 48 with m u / d = 
0.2m s and a « 0.15fm generated by the MILC collaboration has been used. As a basis one qq 
and one meson-meson interpolator, picking the goldstone pion/kaon K5/K5 is used. Extracting two 
energy levels and assuming a Breit-Wigner shape the resonance parameters in Tableware extracted. 



1=\ kK 


7 = Kit 


g^ = 4.54(76)GeV 
M R = 0.779(27)a 


g ann = 2.69(44)GeV 
M R = 0.691 (37)a 



Table 5: Resonance parameters attributed to the k and a resonances [|5^, |56p. For cautionary remarks please 
refer to the text. 



It should however be pointed out that the qq interpolators inevitably couple to all staggered 



taste combinations Q5j> |58|]. Therefore, the variational analysis may render other taste combinations 



Kbftb an d ftbitb as excited states. Moreover the experimental data already shown in Figure [13| for 
the K channel makes the assumption of a Breit-Wigner shape questionable. 

3.2.3 A(1232) ^ Nk by the QCDSF collaboration 

The QCDSF collaboration is investigating nucleon-pion scattering in the A(1232) channel 



[p9p using 32 3 x T, 40 3 x T and 48 3 x T lattices with a pion mass m K w 250 MeV. Preliminary 



results for the phase shift are plotted in Figure 17. For the delta-pion-nucleon coupling they obtain 



2 

= 11 ±3 which should be compared to the value A7 ^ K exp w 14.4 extracted from experiment. 

3.2.4 Baryon interactions in a matrix Hamiltonian model 

In Luscher's method finite volume errors should be exponentially suppressed as a function of 
the lattice size L. This assumes that hadrons at the boundary of the box are in the asymptotic region, 
otherwise the error may scale as L _1 . In [EO] this issue is investigated in a matrix Hamiltonian 
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Figure 17: Preliminary results from a simulation of nucleon-pion scattering by the QCDSF collaboration 
[^9|]. The plot shows data in the vicinity of the A(1232) resonance. 
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Figure 18: Left panel: Estimates of the phase shift employing Luscher's formula to the ANn model for 
different lattice extent L; Right panel: The resonance energy from the Luscher method compared to a direct 
calculation. Plots taken from 



model for A Nn. Figure [l8| shows some of the results. In the left panel results from the Luscher 
method using a different spatial extent are plotted. In the right panel the results from the Luscher 
method are compared to the direct calculation. While both methods agree for large volumes, it is 
clear that large finite volume effects are present in this case. For further discussion please refer to 




4. Beyond QCD 

As already demonstrated in the case of toy models, nonperturbative methods for the extraction 
of resonances in a lattice simulation can also be applied to problems beyond QCD. As an example, 
the Luscher method has been used recently in a simulation of the pure Higgs-Yukawa sector of 
the electroweak standard model []6l|]. The model consist of a complex scalar Higgs doublet and 
a mass-degenerate fermion doublet (representing top and bottom quarks) coupled in a chirally 
invariant way. In this case properties of the Higgs resonance can be extracted. 

Figure [l9| shows the results for three different quartic couplings X. While this model is quite 
distinct from the situation in nature, is is very interesting to observe that the width of the Higgs 
stays small even at large coupling. 
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Figure 19: Elastic scattering phases obtained for three different quartic couplings A within a chiral Higgs- 
Yukawa model. Plots taken from [^TJ. 

5. Summary & Outlook 

While the field is progressing fast and the last year has seen first QCD resonance studies in 
channels other than the p meson channel, lattice studies of (QCD) resonances are still in their 
infancy. Part of the reason is that the standard Liischer method is limited to the case of elastic 
scattering. While there has been some progress recently []62|-|66|] most ideas to go beyond the case 
of elastic scattering require a certain degree of modeling. 

To illustrate the severity of this limitation it is instructive to appeal to experiment once more. 
Even in the low-lying meson spectrum, there are several interesting cases where multiple thresholds 
are expected to be important. Two prominent examples would be the ^o(980) which goes both to 
T]7T and KK and the ^i(1270) where the branching ratios to Kp, K ic (872)n 9 Kco and ££(1430) ft 
are all known to be sizable. 

Last but not least there are also many interesting states, especially the so-called X, Y and Z 
states in the charmonium spectrum, where a first principle calculation from lattice QCD is needed, 
but where most interesting states are above multiparticle thresholds. In short, there is much progress 
to be made for a comprehensive description of QCD resonances, and beyond. 
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